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Abstract 

We study the dynamics of dissipative spins for general spin-spin coupling. We investigate the population dynamics and 
relaxation of the purity in the white noise regime, in which exact results are available. Inter alia, we find distinct reduction 
of decoherence and slowdown of purity decay around degeneracy points. We also determine in analytic form the one-phonon 
exchange contribution to decoherence and relaxation in the ohmic quantum noise regime valid down to zero temperature. 
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1. Introduction 

The spin-boson model is a key model since the 80's 
for the quantitative study of decoherence, relaxation 
and energy dissipation [1,2]. In view of the substantial 
progress in fabrication of coupled qubit devices [3,4] 
and major advances in quantum state manipulation 
towards quantum computation [5] , there is growing in- 
terest in the accurate calculation of the dynamics of 
coupled spins for realistic environmental couplings. 

Here we communicate new analytical results for the 
dynamics of two mutually interacting spins each liable 
to independent white noise or to quantum noise forces. 
Previous work [6] is extended in three different direc- 
tions. First, we determine the dynamics of the purity 
in the white noise regime (WNR). Second, we study 
the dynamics near two different degeneracy points and 
find striking reduction of purity decay and decoher- 
ence. Third, we calculate the exact one-phonon con- 
tribution to dephasing and relaxation in the quantum 
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noise regime in analytic form. The resulting expres- 
sions hold down to zero temperature. The major results 
are obtained within the real-time path sum method for 
the 16 states of the two-spin density matrix. The en- 
vironmental couplings are included via the Feynman- 
Vernon method. Here we omit methodical and techni- 
cal aspects. We rather put emphasis on results in ana- 
lytic form and their physical implications 

In section 2 the model of two interacting spins each 
coupled to its own dissipative environment is intro- 
duced. An analytic analysis of the white noise regime, 
inter alia study of purity and degeneracy points, is 
given in section 3, while section 4 deals with the effects 
of weak quantum noise. 



2. Model 

We consider two two-state systems or spins which 
are mutually coupled via Ising, XY, and/or Heisenberg 
coupling. In addition, each of them is coupled to its 
own heat bath. In pseudospin representation, the two- 
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spin-boson Hamiltonian reads (we put R = fee = 1) 

H = Hss + HsR . (1) 
Here, Hss represents the interacting two spins, 

HsS = —5A1 CTj, — iA2ra - ^ 



jCl (Tz — o£2 Tz 



(2) 



In the basis formed by the locahzed states \R> and \L>, 
the parameters 2 ^nd Ai,2 represent the bias ener- 
gies and tunneUng couphngs of the a- and T-spin, and 
Vx,y,z are the interaction parameters. The term Hsn de- 
scribes the spin-reservoir couphngs and the reservoirs, 

-ffsR = - - -TzX2 + J2 Z^'^C,a4,a^C,c« • 

C = l,2 a 

Here, X^t) = E„ Cc,a[&c,a W + bl„{t) ] (C = 1, 2) is a 
collective reservoir mode. All effects of the environment 

arc carried by the power spectrum of the collective bath 
modes. We have 5^ <./(cij) — S^^^,Sq{u!), where 



/oo 
dte"^\X^{t)X^{Q))^ 
-00 

= 7rGf(a;)coth . 
The spectral density of the coupling is [1,2] 



(3) 



kl/<^c 



The second form describes ohmic coupling with high- 
frequency cut-off and dimensionless damping con- 
stant Kf. Instead of the independent baths, one might 
also choose a common bath for the two spins [7] . 

In the sequel, we confine ourselves to the case of Vy- 
and Wj-coupling of the two spins. This case is most in- 
teresting concerning application to coupled Josephson 
junctions. In addition, we disregard the bias terms. 

The Hamiltonian Hss is diagonalized with the uni- 
tary matrix 



f cos (01 ) 

/+(0i) 

cos {(j>2) 



U ■■ 



cos (02) 
/+(02) 
cos (01 ) 

V /-Wi) 



/+(<^l) 

/-(02) 
-/+( 

-/-( 



/+(<^l) 
-/-(02) 
2) /+(02) 

1) -/-(<^l) 



cos (01 ) \ 



/+ 




COS 


(<^2) 


/- 


{<i>2) 


COS 


(02) 


/+ 


(02) 


COS 


(0l) 


/- 


(0l) 



with f±{(f)) = x/1 ± sin0 
arctan [{vy =F i!z)/(Ai ± A2 



and mixing angles 0i,2 
)]. We then get 



H = UHU-'' = -^{az^l)-^{l^Tz), (4) 



with the eigen frequencies 

n = i(n+-Fn_), 5 = i(n+-n_), 
fi± = Vi^i ± ^2)2 + {vy T vz)2 . 



We have the Vieta relations 

nl+if_=2 (A? -h A^ ^ Uy 



2 , 2^ 
+ Vz) 



(5) 



(6) 



il^ (5^ = {vy Vz — Ai A2)^ 



The two-spin density matrix has 16 matrix elements. 
They can be expressed as linear combinations of the 
unit matrix and the following 15 expectation values, 
{ai (g) l)t = {ai)t, (1 <8) Ti)t = {Ti)t, and {at (g) Tj)t = 
{<^iTj}t {i = X, y, z and j = x, y, z). These quantities, 
denoted by Wj{t) {j = 1, - ■ ■ , 15), obey the equations 
of motion Wj^t) = i[Hss, Wj{t)] {j = I,-- - ,15). 
The set of coupled equations are conveniently solved 
in Laplace space A. Throughout we choose the initial 
state {crz)o = 1, {tz)o = 1, and {azTz)o = 1, and all 
other expectations zero. 

The resulting expressions may be written as 

WjiX) = Nj{X)/Dj{X) , i = 1, • • • , 15 , (7) 

where the AT, (A) are different for the individual Wj{X), 
while the Dj{X) fall into the following two categories, 

A2,^(A) = {X^ + n''){X' + S'') , 



Dn±{X) = X{X^ + nl){X'' + . 



(8) 



3. White-noise regime 

3.1. General features and qualitative behavior 

The exact formal solution of the dissipative two-spin 
dynamics has been discussed in Ref. [6]. Explicit ex- 
pressions for the Wj{X) have been given in the white 
noise regime (WNR), in which (3) reduces to 



SciuJ < T) = 2i?c 



where i9 



2nKcT 



(9) 



It was found that the form (9) is expedient in the regime 
^ 0.3 and ii± T -C Wc, which covers not only the 
incoherent regime but also a sizeable domain of the co- 
herent regime. The analysis shows that reservoir modes 
in the range 2ttT < u) < lo^. give rise to an adiabatic 



2 



(Pranck-Condon-typc) renormalization of the tunnel- 
ing coupling, ^ = (27rT/A<;,r)^^« A^_, with 
aJ j.^"- = Af/wf''" [2,6]. In the reminder of this sec- 



tion, we assume that the Ai,2 are the renormalized 
ones. The modes with cu < 2nT lead to dccolicronce 
and relaxation. In the WNR, they are accounted for 
by an appropriate shift of the Laplace variable A in the 
time interval, in which spin C, dwells in an off-diagonal 
state. We have (C = 1,2) 

A ^ Af = A + , and A ^ A12 = A + t9i + i?2 . 

The resulting 15 coupled equations are 



A (cr^) 


lA/\ / \i/ \ 

= 1 — Ai(CTy) - Vx{TxCry) + Vy{axTy) , 


Al {(^y) 


= Al(az) + Vx{(JzTx) — Vz{Tzax) , 


Al {a^) 


= —Vy{crzTy) + Vz{(TyTz) , 


A(r.) 


= 1 — A2{Ty) — Vx{CxTy) + Vy{ayTx) , 


A2 {Ty) 


= A.2{Tz) + VxipxTz) - VzipzTx) , 


A2 {Tx) 


= -Vy{cTyTz) + Vz{cTzTy) , 


A12 {UxTy} 


= A2{axTz) + Vx{Tz) — Vy{(Jz) , 


A12 (lyyTx) 


= Ai{Tx(Jz) + Vx{(7z) — Vy{Tz) , 


M{iyzTy) 


= A2{azrz) — Ai{ayTy) +Vy{ax) — Vz{tx) , 


XliOyTz) 


= Ai{azrz) — A2{ayTy) +Vy{Tx) — Vz{(Jx) , 


Al {UxTz) 


= —A2{axTy) — Vx{ry) + Vz{(Ty) , 


A2 {(JzTx) 


= —Al{rxay) — Vx{ay) + Vz{Ty) , 


A12 {OyTy) 


= Ai{azTy) + A2{ayTz) , 


A {azTz) 


= 1 — Ai{ayTz) — A2{azTy) , 


A12 {OxTx) 


= 0. 



The solutions for all Wj{X) are again in the form (7). 
There arise only three different denominators, namely 

Di(A) = (AAi + A?)(AiAi2 + Ai) + vl{Xl + v^) 

+ AAi2i'2 — 2vyVzAiA2 , 

i:'2(A) = (AA2 + Ai)(A2Ai2 + Al) + vliXl + ) 
+ AAi2V^ - 2vj,v^AiA2 , 

DsiX) = AAi2[A?Ai + (A? + Xl) + 2v^AiA2] 
+ (t-^A^ + v^.Al + A^AiA2)(AiAi2 + AA2) 



+ {vlAi + vlAl + A?AiA2)(A2Ai2 + AAi) 
+ {vl - vlfXX^2 + {Al - Alf\iX2 
+ 2wj;i;;,AiA2(A + Ai2)(Ai + A2) . 



For instance, we obtain 

Al (Ai + AiAi2) -|-w?Ai2 



(c^.(A)> = 
{rz{X)) 



DiiX) 

A2 {Al + A2A12) + i;^Ai2 



152(A) 

{ayTy{X)} = [AiA2{vlWz + A1A2) 

, / A 2 , A 2n 1 Al + A2 

+ VyVz[Ai+A2)\ — 



(10) 



D,{X) 



For vanishing bath coupling, both Di (A) and D2 (A) go 
to Dq^s{X), and D3{X) — » A reduces to Dn^ (A). 

The dynamics of the expectations Aj{t) is mainly 
determined by the zeros Ai ofZ)fe(A) (fc = 1, 2, 3). They 
appear in complex conjugate or real pairs. We get 



(11) 



where n is either 4 or 6. The behaviors of the four Xj of 
(cz(A)) (and of {tz{X))) and the six A^ of {ayTy{X)), and 
the respective amplitudes Bj, are quite diversified. In 
Fig. 1 we show plots of the real parts (rates) and imag- 
inary parts (oscillation frequencies) of the four Xj of 
{(Tz{X)} as functions of i9 for a particular set of param- 
eters (identical spins, Ai,2 = A, i?i,2 = 1?). In the cou- 
pling range Vz < v* there are three crossover tempera- 
tures 1?! , and -dfi at which the discriminant of Di (A) 
is zero. For Vz > there is only a single crossover 
temperature 1^3. The critical coupling strength is 

v: = l[{2A'+vl)'/'-Vy]. (12) 

A plot of the crossover temperatures is shown in Fig. 2. 

The particular case Ai = A2 and Vz = —Vy is a 
degeneracy point, Q = 5 (cf. subsection 3.4). In this 
case, the crossover curves i^Kvz) and ^2{vz) coincide. 

In the regime -d < '&2 the dynamics is coherent, de- 
scribed by a superposition of two damped oscillations 
with amplitudes of comparable size. For 1? < the 
two oscillations have different frequency, Q > S, but 
the same damping rate. In the range '&{ < -& < 1&2: they 
have the same frequency, but different damping rates. 

In the range 1? > 'i?2 , the dynamics is incoherent with 
4 different relaxation rates. The two smallest rates have 
sizeable amplitudes and thus dominate the relaxation 
dynamics. In the so-called Kondo regime i9 > 1^3, three 
of the four Aj-contributions have negligibly small am- 
plitudes. The only relevant contribution is that with 
the smallest rate. The Kondo characteristics is that 
this rate decreases with increasing temperature, 7k = 
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Fig. 1. (<Tz): Plots of ReAj and ImAj against i? = i?i^2- 
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Fig. 2. (o"z}- Crossover temperatures against Vz for 
Vy = 0.1 and identical spins, Ai,2 = 1, and i9i.2 = i9. 

A^/i9 (cf. Fig. 1). This counter-intuitive feature is al- 
ready well-known from the ohmic single-spin model [2] . 

3.2. Low temperature WN regime 

The WN regime has a low temperature bound 
roughly given by T ~ Above this bound and be- 
low the first crossover temperature, Q,± < Ti , the 
real parts of the {Xj} vary linearly with temperature. 
In this regime, systematic low-temperature expansion 
of the zeros of Dj{X) {j — 1, 2, 3) is straightforward. 
The results are as follows: 

{o'z)t'- There is a superposition of two damped os- 
cillations, Ai_2 = iif^ ~ 7n ^i^d A3 4 = ±i(5 — 7^. 
The frequencies S7 and 5 are close to their bare val- 



ues given in Eq. (5) near T — i}±. As T is increased, 
they approach each other and coincide at T = Tj* . 
The respective damping rates 7q and 7^ in the regime 
n± <T < Ti read 



1 , n'-Aj-vl ^ Q?-Al-vl ^ 



1 ^ lil+vl~5\ Al+vl~ 5^ 

'ya = + ^^ + — 

2 ' 2(^2-52) ' 2(02-52) 



(13) 



The amplitudes of the oscillatory contributions (in ze- 
roth order in i9i,2) rea 



Mutual exchanges Ai «-» A2 and i9i ^ i?2 yield the 
corresponding rates and amplitudes for (r^)*. 

{o'zTz)t'- According to the zeros of D3(A), Ai,2 ~ 
±ifl+ -7n^, A3,4 = ±iQ- -7r!_, Ag^g = -75,6, there 
are two damped oscillatory and two relaxation contri- 
butions. The frequencies fl± are close to their bare 
values given in Eq. (5) near T — Q± , and they coincide 
at the first crossover temperature Tf. The damping 
rates 7^^ and amplitudes Bqj_ of the oscillations are 



7n+ = 7fi_ = 5(^1+^2), 



Br 



(Ai A2 
4^2 



Bn_ 



(Ai-A2)^ (14) 

4nl 



The relaxation rates 75_g are determined by a quadratic 
equation, which is obtained by truncation of 1)3 (A). 
The resulting expressions for the relaxation rates and 
associated amplitudes are 



75 



76 



- A? - „ fi^ - Aj 



02-52 

Al+vl-S^ 

02-52 

{VyQ, + Uz5)^ 

(02-52)2 



i9i + 



02-52 



Al + vl- 52 



(15) 



02-52 



(02 - 52)2 



3.2.1. The limit A2 — » and Vy ^ Q 



In the limit A2 — + and 1"^— +0, the transition frequen- 



cies 0± become degenerate, and hence O = \J A\ + ii| 
and 5 = 0. Accordingly, the expression (13) for 7^- is 



The amplitudes ,5 and are the residues of each 

of the corresponding two complex poles. 
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not valid anymore. To cope with this limiting case, we 
must determine A3, 4 from a quadratic equation, which 
is found from Di (A) = by reduction. The respective 
complex eigenvalues for slight detuning are found as 



A = ±i\ 5^ 



4^n2 

From this form we see that the two complex conjugate 
eigenvalues turn into two real eigenvalues, when S is 
sufficiently small. At 5 = 0, the rate expressions are 

7+ = ??i + i?2 , and 7_ = 7^ := i?i . (16) 

In addition, the analysis shows that the residuum as- 
sociated with the pole at A = — 7_|_ is zero, while the 
other yields the amplitude = ivf/Jl^. 

At this point, we remark that, in the limit A2 — » 0, 
the coupling Vz takes the role of a biasing energy for 
spin a. Thus, the dissipative two-spin model reduces 
to the standard biased spin-boson model. The rate 7r 
is just the relaxation rate of this model. Furthermore, 
the rate 7^ in Eq. (13) reduces to the form 



7n = H + {vl/n^)i). 



(17) 



This expression coincides indeed with the decoherence 
rate of the biased spin-boson model in the WNR [2] . 

3.3. Purity 

For a system described by the density matrix p(t), 
the purity P{t) := Trp^{t) tells us whether the system 
is in a pure state or in a mixture. For a pure state, 
there is P = 1 while for a fully mixed state P — 1/N. 
Here N ist the number of the system's accessible states. 
In the low temperature WN regime discusssed in the 
preceeding subsection, the purity P{t) is found as (the 
index a, r refers to the respective spin) 



P(t) = i + i(l-a 



'* + i(l + a)e-'^".<' 



+ 1(1 - Cr)e-^^^.-* + 1(1 + CO e-'^n 



+ i(l-Co_)e 



-27n 



+ i(l-C7oJe 



+ 



4 (J^2_J2)2 

with the amplitudes 



"^4 (02 _ 52)2 



A? 2 - Ail + vl 



^2 -52 



Cn = 



This function smoothly drops on the time-scale given 
by the system's damping and relaxation rates from the 



initial value P{Q) = 1 to the fully mixed thermal equi- 
librium state, P{t — > 00) = i. Observe that alldephas- 
ing and relaxation rates relevant to the decay of the 
expectation values Wj{t) (j = 1, ••• , 15) contribute 
to the decay of the purity. 

3.4. Decoherence dip near degeneracy points 

Of particular interest are degeneracy points of the 
two-spin system. There are two different cases: 

Ai = A2 , Vy = —Vz =^ Q, = 5 (case 7) , 



VyVz = A1A2 



(case 77) . 



For comparison, we also study the nondegenerate point 
conjugate to case 7 



Ai 



n,<5 = Ait; (case 7*). 



Consider first ((Jz)t in case 7. We find from Eq. (13) 
upon taking the limits Ai,2 = lim^^o A ± i?; and 
Vy^z = limK-.o -iiv+^K the rate expressions 

ln,8 = \&i + \&2 T + 02) T - ^2) , 



where Q. = ^/A^ + ^P. These forms reduce for equal 
bath coupling ■& = i?i,2 to 



1? T ^1 



(18) 



Thus, the one rate is smaller and the other larger than 
lO. The amplitudes associated with (18) are found as 



Bn,5 



lA^ A + . 
4 V !^ 



(19) 



Hence in the regime w <C A, the amplitude of the 
smaller rate is maximal, while that of the larger rate is 
negligibly small. 

The results (18) and (19) may be compared with 
those of the nondegenerate case 7* . They are 

70 = 75 = I i5i -h 3 i?2 , and Sn = = i . 

The decline of decoherence at the degeneracy point 
7 compared to point 7* is clearly visible in Fig. 3. The 
decoherence minimum follows from competition of the 
two equally preferred ground states. Due to the Vy- and 
Uj-couplings, the system could relax either to parallel 
or to antiparallel alignment in y- or z-direction. Hence 
the respective second spin-spin coupling gives rise to 
partial suppression of decoherence. Reduction of deco- 
herence at point 7 may be looked upon as a new type of 
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Fig. 3. Plots of {cr^)t against t for identical spins at the de- 
generacy point / (full curve) {vy = —Vz) and at the nonde- 
generate point I* , Vy = (dashed curve). The parameters 



frustration of decoherence [8] . Here, the phenomenon 
is due to the non- commutative spin-spin couplings. 
For {azTz)t the picture is similar. In case /, we have 
— 2Q, Q- — 0, and the rates and amplitudes read 



^7n_ = |(i?i+^2) 



75,6 = T; =F T^i^l - ^2) =F + 1^2) : 



2 ' 20. 
^ A2 1 



R - R - " 

'+ 4^2 ' 4 ' " ^6 " 4^2 

Thus we find for equal bath couplings i9i,2 ~ 1? 



75,6 



(20) 



In contrast, in the non-degenerate case /* we have 
— 2A, f2_ = 2u, and the rates and amplitudes are 

7n+ = 7n_ = + i?2) , 15,6 = U^i + ^2) , 
Bn^ = I , Bq = , = = i . 

Thus we have decline of relaxation of {azTz)t at 
the degeneracy point /. In Fig. 4 we show a plot of 
— -fTTTv — , which is a form of an effective transition 

P(t) at ' 

rate from a pure to the fully mixed state. At the de- 
generacy point C, = —1 (case 7), there is a distinct 
slowdown of the extinction of the pure intial state. 

Case 77 is another set of parameters for which the 
spectrum is degenerate, r2+ = 57 _ [9]. An expedient 
parametrization for identical spins, Ai_2 = A, is 

2 l,-2 2\ 2 1/2, 2\ 

V„ = i;(v —V_\. V, = iv +V_^]. 

2 2^2 . r^-^ . (21) 

Then we have 



y^v* - 4A4 > 0. 



= 2A^ + 



S = . 

The undamped dynamics of {(Jz)t is 
^2 I ^2 

(a,)t = cos(nt) + ~ [1 - cos(nt)] 



(22) 



Fig. 4. The quantity 



is plotted against 



P(t) dt 

C = I'y/^z for identical spins. There is a distinct lowering of 
the effective damping rate at the degeneracy point f = — 1 
(case I). The parameters are A = 1, K=0.01 and T = 3. 

as follows from Eq. (10). At the degeneracy point 77, 
the rate expressions (13) would yield 

ln,s = \{M^+^2)±^{^^-^2). 
For identical bath couplings, these would reduce to 

1^ = 16 = ^- (23) 

Now, as we have already argued in subsection 3.2.1, 
the expressions for 7^ are not correct near to and at 
the degeneracy point. In fact, then the poles of ((Tz(A)) 
at A = ± i 5 — 7^ are determined by a quadratic equa- 
tion in A. As one approaches the degeneracy point, the 
complex conjugate roots turn into two real ones. We 
find in the limit 5-^0 for i9i,2 = 



T i <5 + 7i 



7± 



2^2 



Furthermore, the respective amplitudes are 

^" = 4(^-?p)' ^-. = 4(i + ?7^)^2n- 

Now, since the rate 7_ is smaller than and the re- 
spective amplitude is nonzero, we have again, now at 
point 77, reduction of the decay of the purity P{t). 

Similar behavior occurs also in {pzTz)t- While 7^^ 
and 7q_ are as in Eq. (14), we now have for i9i,2 = 1? 



75 



76 



and the amplitudes read 



Bn_ = 0, 



B^ = 

'5,6 



2r22 



The decrease of the rate 7g for > and of 75 for 
> leads again to a slowdown of the decay of P{t) 
around the degeneracy point 77. We have verified nu- 
merically that in the quantum noise regime T <^ A 



considered below, in addition to the rate 75 or 75 , also 
the dephasing rates 7^^^ and 7^^ exhibit a pronounced 
minimum at the degeneracy point //. As a result, be- 
sides the indentation in the purity characteristics (cf. 
Fig. 4), also the dephasing of the two-spin dynamics is 
considerably slowed down at the degeneracy point II. 
A numerical analysis of this phenomenon at the degen- 
eracy point // is reported in Ref. [9] . 



4. Quantum noise regime 

Consider next the extension of the analysis to the col- 
ored quantum noise regime (QNR) relevant atT^ f2± . 
Since at low T quantum noise prevails, Eq. (9) is not 
valid anymore. Rather we have to revert to the expres- 
sion (3) . In the ohmic case, we have 



2ttK( uj coth 



\2T) 



(24) 



We have studied the effect of the one-phonon exchange 
contribution to the dynamics of the two-spin model 
using both the perturbative Redfield approach [10] and 
the self-energy method within the path sum method 
[2]. The latter amounts to systematic calculation of 
the self-energy to linear order in the bath correlations, 
Re{X^{t)X^{Q))T = Qc(t) (C = 1, 2). We have 



Qdt) 



i r 

1" Jo 



[ 1 — cos(a;t) ] . 



(25) 



In the standard notion of sojourns and blips [1], the 
one-phonon self-energy, say Si (A), receives contribu- 
tions from the intra-blip correlation and from the four 
intor-blip correlations induced by bath 1 between a pair 
of blips of the cr-spin. The inter-blip correlations vanish 
in the WNR. In the usual charge picture, the former 
correlation is a charge-charge interaction and the latter 
correlation corresponds to a dipole-dipole interaction. 

In the time intervall between the correlated blips, 
the cr-spin may perform any number of uncorrelatcd 
jumps between its two blip and two sojourn states. 
In addition, we must take into account all transitions 
which spin t can make during the dwell time of spin a 
in sojourn and blip states which are spanned by bath 
correlations. The succession of flips of the a- and t- 
spin is dictated by the Hamiltonian (1) with (2). Fol- 
lowing the lines expounded for the single spin-boson 
model [2], it is straightforward, but tedious, to calcu- 



late the self-energy Ei (A). Interchange of the two spins 

and reservoirs then yields S2(A). 

The self-energies Ei,2(A) lead to shifts of the poles of 
the Wj(X). It is advantageous to measure the resulting 
shifts in terms of generalized scaled temperatures 0i,2. 
These depend on the power spectra (3) and are nor- 
malized such that they reduce to the previously intro- 
duced scaled temperatures i?i,2, Eq. (9), in the white- 
noise limit. For lack of space, we now put Vy =Q. 

{cz}t- The damping rates of the two oscillations with 
frequencies Q and S are found to read 



7r! = 



7« = 



where 



2^2 - A? 



2 (n^ _ 52) 
2 (^2 - 52) 



S2) ^2 



(n) 



2(n2 

2 (02 _ 52) 



(26) 



- 2 



TT 2 (^2 - Al)Si{S) + (A? - 52)5,(J7) 



2n^ - 52 _ Af 

TT {Q^ - Ai)Si{5) + 2{Al - 5^)Si{Q) 



fi2 + A2 - 252 



ef' = i7r52(fi), 



e 



i7r52(<5). (27) 



The amplitudes asociated with the complex frequencies 

A = =Fi n — 7q and A = =Fi 5 — 7^ are 



Bn = 



Q?-Al-vl „ vl + Al-5^ 

-fc><5 = 



2(^2-52) ' 2(^2-^2) 

These onc-plionon rate expressions hold in the QNR 
down to T = and they smoothly map on the WNR 
results (13) at elevated temperatures. 

In the corresponding expressions for (T^)t, the in- 
dices 1 and 2 are interchanged. 

Following the lines expounded in subsection 3.2.1, 
we may also consider the limit A2 0. In this limit, 
the characteristics of the pole trajectories is as in sub- 
section 3.2.1. The resulting forms for 7^ and 7[j are 
those of the biased spin-boson model in the one-phonon 
QNR [2], 

,,2 



7^1 + 1^51(0). (28) 



Consider next the limit — > 00 for the symmetric 
case Ai_2 = A and A'1,2 = K- For large coupling, the 
two spins are locked together and behave like a single 
spin with oscillation frequency 5 = A^/d^. Since the 
amplitude Bn becomes neglibly small, the dynamics is 
(cT^)t = cos(5i) exp(— 7jt) with the dephasing rate 



7s 



■n K 5 coth 



(2t) ' 



(29) 
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as follows from Eq. (26) . Since the effective spin is un- 
biased there is no relaxation term. In the WNR limit, 
the rate (29) reduces to 7j = i?, which is twice the de- 
phasing rate of the spin-boson model, Eq. (28) , in this 
limit. The additional factor two is because the effective 
spin is coupled to two identical reservoirs. 

The corresponding expressions for {Tz)t follow from 
these forms by interchange of the indices 1 and 2. 

{pzTz)t'- As temperature is lowered from the WNR to 
the QNR, the damping rates of the oscillations with fre- 
quencies f2± change from '^^j^ = |(^?i + 1^2), Eq. (14), 
to the one phonon expression 



7r!4 



0' 



(n±) 



2 



e 



(n±) 



^2 

2 n2 



2 02 _ 52 



&(n). 



5. Summary 

We have studied the dynamics of a spin or qubit cou- 
pled to another spin. The latter could be another qubit, 
a bistable impurity, or a measuring device. Wc have 
given the dynamical equations in the WNR for general 
spin-spin coupling and we have discussed the rich fea- 
tures of the coupled dynamics. Analytic expressions for 
dephasing and relaxation rates and for the decay of the 
purity have been given in the one-phonon WNR limit. 
Furthermore, the corresponding generalization to the 
quantum noise regime, which is based on a systematic 
calculation of the self-energy, has been presented. Our 
results smoothly match with those of the perturbative 
Redfield approach. 

Financial support by the DFG through SFB/TR 21 
is gratefully acknowledged. 



As regards the relaxation rates 75 and 7g of {azTz)t, 
the situation is more subtle, because they are deter- 
mined by a quadratic equation in A which involves the 
self-energy Ei,2(A) in linear and second order in ii'1,2. 
The calculation is most easily performed within the 
Redfield approach. The resulting rate expressions are 



76 



76 = 



n2 

<2 



^1 c,(0) , 

"75" ^1 + Tt; 



^2 ^(0) 



^1 - ^ q(o) , ^2 

n2 - 52 ^1 + 



(30) 



02 - J2 "^2 : 



and the amplitudes read 

^T5 - (^2 _ 52)2 ' -"Te - (J^2 _ 52)2 ■ ^^^) 

The functions @^^^2 depend on the power spectra at the 
transition frequencies Q. and 5. With the abbreviaton 
Af ± A| = Aj- , we find the explicit form 



©1°2 = J{ft,2(n)+ft.2(5) 



{Q? - 52)A2 
2«'Ai,i 



[SiA^) - 5-1,2(5)] 

[&.i(n)-5'2.i(5)]} 



(n2-,52)A?. 

These expressions hold under assumption Ai l^i. 
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